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We analyze numerically a two-dimensional \<f> 4 theory showing that in the limit of a strong 
coupling A — > oo just the homogeneous solutions for time evolution are relevant in agreement with 
the duality principle in perturbation theory as presented in [M.Frasca, Phys. Rev. A 58, 3439 
(1998)], being negligible the contribution of the spatial varying parts of the dynamical equations. 
A consequence is that the Green function method works for this non-linear problem in the large 
coupling limit as in a linear theory. A numerical proof is given for this. With these results at hand, 
we built a strongly coupled quantum field theory for a Xcj) 4 interacting field computing the first 
order correction to the generating functional. Mass spectrum of the theory is obtained turning out 

' to be that of a harmonic oscillator with no dependence on the dimensionality of spacetime. The 

agreement with the Lehmann-Kallen representation of the perturbation series is then shown at the 

, first order. 

(N ; 

A lot of problems in physics have such a difficult equations to solve that the most natural approach is a numerical 
one. Weak perturbation theory generally proves to be insufficient to extract all the physics. A well-known case is 
given by quantum chromodynamics that due to the strength of the coupling constant at low energies, makes useless 
known perturbation techniques demanding the need for numerical solutions. 

In the seventies and eighties of the last century a significant attempt to build a perturbation theory for a strongly 
' interacting quantum field theory was proposed plStliEllHISHIl]- ^ n * ms approach it was stipulated that the 
perturbation to be considered is the free part of the Lagrangian. Notwithstanding this approach is still studied today 
qq ■ [9j no fruitful results have been obtained so far due to the strongly singular perturbation series that is obtained in 
^sO [ this way. Rather, the rationale behind this method is really smart as one recognize that just interchanging the two 
■ parts of the Lagrangian one gets different perturbation series. 

This duality in perturbation theory is a general mathematical property of differential equations as was shown in 
Ref . [ToL HH] . What makes duality interesting is the general property of the leading order that, while in the weak 
perturbation case is just a free linear theory whose solution is generally known, for the dual series that holds in the 
limit of a strongly coupling, that is a coupling going to infinit y, o ne can prove a theorem showing that the adiabatic 
approximation applies. We also pointed out in recent works [P2l Il3| that in field theory and general relativity the 
dual perturbation series at the leading order produces a rather interesting result: in a strongly coupled field theory 
the leading order is ruled by a homogeneous equation, that is, the spatial variation of the field in the equations of the 
theory becomes negligible. In general relativi ty t his give s precious informations on the space-time near a singularity 
where the above behavior was conjectured in [bH Il5l Tl6| and numerically shown in |l7j . 

In this paper we have two different aims. Firstly we intend to prove that the numerically observed behavior in 
general relativity also holds for a A</> 4 theory, that is, a homogeneous equation rules the leading order of a strongly 
interacting scalar field. Then, after numerically proving that in a strongly coupled field theory the Green function 
can be used in the same way as done in a weak field theory, a quantum field theory is obtained. 

We apply the duality principle in perturbation theory as devised in 0, 0, fl^L Il3j assuming a Hamiltonian for the 
field (here and in the following we take h = c = 1) 
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H = I d°- x x 
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being D the dimension, the mass and A the coupling. For our aims we take (j,q = 1 and a single component scalar 
field. This Hamiltonian gives the following Hamilton equations 

d t 4> = vr (2) 
<9 t 7T = d 2 x (j) - <j> - XV \<P) 

apex meaning derivation with respect to <f>. From eqs.© we recognize two perturbation terms d^cj) — <f> and V'((f>) and 
one may ask what is the relation between the weak perturbation series for the latter term with the one having the 
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term 9 2 — as a perturbation. Indeed, by exchanging d 2 — «-> V'(4>) for perturbation the following equations can 
be obtained 

<9 T 0o = tto (3) 

5 T 0i = 7Tl 

a r7 r = -f'(^o) 

9 T 7n = -F"(0 o )0i + - 0o 



This is a non trivial set of equations that can be recovered if we take 
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So 



our interchange of the perturbations produced a dual series that holds in the limit A 



produced a dual series tnat nolds m tnc limit A — > oo as expected by the 
duality principle in perturbation theory |Tol . . The most important result we have obtained is that we get at the 
leading order a homogeneous equation, that is, a self-interacting scalar field with a coupling constant going to infinity 
is ruled by a homogeneous equation. This result is relevant as settles the physical meaning of homogeneous solutions 
for a given field theory. 

Now, let us specialize the above analysis to a A0 4 theory. When A — > oo we have to solve the leading order equation 



-\<t>o that has the following solution by Jacobi elliptic functions 0] 



0o = (2Cx)< 



(VXi + C 2 ), 



(5) 



being sn the snoidal Jacobi elliptic function, C\ and C2 two integration constants that can depend on spatial variables. 
So, this analytical solution has to coincide with the numerical solution of the equation D0 + + \<fi 3 = with A 
very large, after the proper boundary conditions are set. Another interesting problem is to see how farther can be 
considered to hold the approximation 



(j>o{x,t) 



G(t-t')j(x,t')dt' 



(6) 



as a solution of the equation D<p + (f> + A0 3 = j in the limit of A very large and G(t — t') the Green function given by 
the equation d?G(t) + XG 3 (t) = S(t) that is 



G + (t)=6(t)( - 



sn 



The time reversed solution 



GL(t) = -e(-t) ( - 



sn 
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(8) 



also holds. It is not difficult to verify that G-(t) = G+(—t). 

The first numerical analysis we worked out is to verify that indeed, when A is very large, a good first approximation 
is given by the leading order solution (Jsj- In order to check this we consider the equation 0<j) + (f> + A0 3 = for D = 2, 
A = 10 4 and take 0(0, t) = 0, 0(1, t) = 0, 9 t 0(a;,O) = and 0(z,O) = x 2 - x. The solution is given in fig. □ The 
analytical solution in this case can be easily computed by eq.JSJ) giving 



(x 2 — x)sn 



(x - x)\ -t + x ,i 



(9) 



FIG. 1: Numerical solution for U<j> + + \<j> 3 = with A = 10 4 . 



z(x.t) 




FIG. 2: Analytical solution of dt4> + \<f> 3 = with A = 10 4 as given in eq.©. 



being xq = cn _1 (0, i) as to have sn(xo,i) = 1. This solution is presented in fig. [3 and the comparison with the 
numerical result is quite satisfactory. Homogeneous solutions drive, in a first approximation, strongly self-interacting 
scalar fields. 

For the next step we have studied the D — 2 equation □</>+</>+ A </> 3 = — sm(27r(x+i)) with the same value for A with 
boundary conditions <j)(0, t) = 0, </>(l, t) = 0, dt4>{x, 0) = and <p(x, 0) = 0. The numerical solution is given in fig. [3]TIic 
analytical solution can be easily computed with the Green function of eq. Q giving 4> w - Jo G+(t-t') sm{2ir(x+t'))dt' 
and the result is given in fig0] and again is quite satisfactory. These results support the other conclusion that the 
Green function method is still useful in a regime of largely coupled scalar fields. 

There is an exterminate literature for quantum field theory (see e.g. 0,|2(],E] for scalar fields). As a convention 
we will use boldface for D — 1 dimensional vectors. Spacetime signature is (+,—,—,—). We start with the standard 
path integral formulation for the generating functional as Z\j] = j[d<j,] e { l $ dD 4^4>) 2 ^^(v^ 2 ~y 2 ~\v(4>)+ 3 <p\} that 
we rewrite as 

Z[j] = J [d(f\e{^ d D x[±(d t 4>) 2 -\v(4>)+j4>}} e {-ifd D x[i ! (v<f ) ) 2 +} s 4> 2 }} ^ 

separating the leading term from the perturbation in agreement with our discussion above. Using our conclusions 
about Green function obtained above one can write down the generating functional, without the perturbation, by the 
Gaussian approximation 



z o[j] = cxp 



% - J d D x 1 d D x 2 j(x 1 )A(x 1 



(11) 



FIG. 3: Numerical solution for □<£ + cf> + \<j> A = - sin(27r(x + t)) with A = 10 4 . 




FIG. 4: Analytical solution with the Green function of eq.0 and A = 10 4 and forcing function j = — sin(27r(a; + 1)). 



from which one can get the Wick theorem. It is easy to verify that sjJyTM^) 



-iA(y 2 - yi) having set 



A(x 2 - Xl ) = S u -\x 2 - x x )[G + {h - h) + G_(t 2 - h)] = A(a?i - x 2 ). 



(12) 



In order to make all the argument self-consistent we derive the generating functional (|llfl from ea. Hl(J|) . The 
existence of the the leading order functional will rely in the end on the existence of the semiclassical approximation 
for the path integral 

Z \j] - / [d4>] e { i f dD *[H a <+'?- xv M +i *]}. (13) 

This can be seen in the following way. Let us apply the rescaling of time r = \f\t. One has 

Zo[j] = J [#]e{ iVX / d "4l(9^) 2 -^)]} e ^/^^ (14) 

that shows that the limit A — > oo corresponds to the semiclassical limit. That is, the system tends to recover a 
classical behavior in the strong coupling limit and all the results obtained above for this case apply. So, we reinsert 
the original time variable t and take 



(15) 
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being 5<j) a small deviation from the classical solution <fi c that satisfies the equation 

4> c + \V'{<p c ) = j. 

Inserting ea. l|15|l into the functional integral i|13[) one has, using ea. i|lt)|) . 

Z \j]=eiI dD ^F[<j> c }. 

being 

F[<j) c ] = e - iX S d D x[v^ a )-^v'(4>.)] [ [d5(j)]e-^ d ° x ^ ^H+h^V'i^H) 2 



(16) 



(17) 



(18) 

We now apply the property that we have proved for the solution of ea. Q16[) . that is, the Green function method still 
applies as in eq.©. This gives back the Gaussian functional (|ll(l taking into account that, in the limit of interest 
A — > oo, after the substitution of the Green function and (JHJ, F[<f> c ] ~ 1. 

A short digression on the Feynman propagator H12[) is needed. Indeed, it is well known that [l8| 



sn(w, i) 



2tt ^ (_ 1 yi e -(n+i) 7r 



_|_ e -(2n+l)7T 



sin 



(2n + 1) 



2K(i) 



(19) 



being K(i) = J Q 2 



v / l+sin( 



1.3111028777 a constant. This means that a Fourier transform gives 



AM = £ 



Br, 



n=0 



u) 2 — ujj, + it 



(20) 



being B n — (2n + 1) ^7^ i^ e -(2>»+i)^ — an d the mass spectrum of the theory in the limit A — > oo is given by 

i 

w n = (n + |) -j^y (-|) 4 that we can recognize as those of a harmonic oscillator. A mass gap computed for n = 

is given by 5s — jRJT) il) 4 ■ This result does not depend on the dimension D but could depend on the number of 
components of the scalar field that we have not considered here. 

It is straightforward to write down the full generating functional to work out perturbation theory. One has 



Z[j] = exp 



'Sj(yi) Sj(y 2 
that, by expanding the first exponential, gives 

Z\j] = l [ l-^Jd D y 1 d D y 2 (~V 2 + l)5 D (y 1 -y 2 )A(y 1 -y 2 ) 

- \\ d D yi d D y 2 (-V 2 + l)S D ( yi - y 2 )I{yi)I(y 2 ) + . . J 



Zo[j] 



(21) 



(22) 



Zo[j] 



being I(z) = J d D xiA(z — xi)j(xi). We realize straightforwardly that there seems to be a divergence as also happens 
for weak perturbation theory. In order to make the computation physically clear we pass to momentum space by a 
Fourier transform as f(k) = j d D xf(x)e lkx and one has straightforwardly 



z o\j] = exp 



d D k - 



(27T) 



D 



j(k)A(k)j(-k) 



So the first integral becomes 



d D yi d D y 2 (~V 2 + l)5 D ( yi - y 2 )A( yi - y 2 ) = / d D kd D h(k 2 + l)A(k)S D (k + kx)6 D (k + h) 



(23) 



(24) 



and we can dispose of the product of Dirac distributions by substituting one of them with the D-dimensional volume 
Vd divided by {2n) D reducing it to Vd J A^jirh (k 2 + 1) / |^A(w) where we have explicitly given the dependence on 
uj to make clear that this integral seems to diverge and a cut-off in k has to be introduced. But we notice that the 
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last integral is nothing else than A(0) = and so, we take this renormalization constant to be zero. So, finally one 
has 



m = 



d°k 



l ~2j (2^ (k + 1 ll(k)Mk)j(-k)A(-k) 



Zo[j] (25) 



that is the result we aimed to. We notice that to recover the proper ordering in A one has to turn back to space and 
time variables and one can see that we are at order A - ' having the product of two Green functions. We have an 
expansion that holds in the strong coupling limit as promised. We see that this series recover the proper dependence 
on k in the propagator in agreement with the Lehmann-Kallen representation [19l |2fJ . This completes the proof of 
existence of a strongly coupled quantum field theory for a A</> 4 model. 

Recently it was shown by Kleinert as very fine results for critical exponents can be obtained with the variational 
method [23, |U |^ but no hint is given on the structure of the solution of the field equations. Here we have built a 
successful approach showing a possible way to find solutions to non-linear quantum field theories in the strong coupling 
limit. We were also able to show that a homogeneous equation rules the dynamics and Green function methods can be 
successfully applied in the strong coupling limit. All this has been supported by numerical results. So, this approach 
can open up the way to exploit analytical solutions where, presently, just heavy numerical work can be accomplished. 
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